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Abstract: The status of lattice determinations of quarkonia and hybrid meson spectra is presented 
and compared with experiment. Both quenched and unquenched results are discussed. Hybrid meson 
decays are considered. 



THIS PAPER discusses the consequences from 
first principles of QCD for the heavy quark- 
antiquark bound states: quarkonia and hybrid 
mesons. This is indeed the area of hadron spec- 
troscopy where the underlying QCD structure is 
most readily discerned. The heavy quarks act 
as sources for the colour fields which provide the 
binding. It was realised early on ||] that a com- 
bination of a short-range gluon exchange compo- 
nent and a long-range confining force is sufficient 
to give a good qualitative description of experi- 
mental data on quarkonia spectra. This can be 
made more precise by explicit lattice QCD eval- 
uations. 



Hybrid mesons are defined as those in which 
the gluonic component is non-trivial. The easi- 
est way to ensure this is to require a spin-exotic 
state (where the J^*^ value cannot be attained 
from QQ alone). For heavy quarks at separation 
ii, the potential energy of these gluonic excita- 
tions can be established directly by lattice calcu- 
lations Q. This enables evaluation of the mass 
and properties of such hybrid mesons. 



Here we summarise the various lattice ap- 
proaches to heavy quark bound states and present 
the latest results. The most thorough lattice 
studies have been conducted in the quenched ap- 
proximation. We discuss, however, evidence for 
differences between these lattice calculations in 
the quenched approximation and those with Nf = 
2 flavours of light sea quarks in the vacuum. 



1. Lattice QCD and heavy quarks 

The simplest way to treat a heavy quark on the 
lattice is to approximate it as static. The heavy 
quark propagators are then trivial to evaluate 
since they are products of time-directed gauge 
links. This enables potentials between such static 
colour sources at separation R to be defined and 
the resulting spectrum of quarkonia can then be 
evaluated exactly from these potentials using the 
Schrodinger equation in the Born-Oppenheimer 
or adiabatic approximation. One advantage of 
this approach is that the continuum limit (lattice 
spacing a — s- 0) can be readily taken. Moreover, 
lattice results at small separation in terms of the 
lattice spacing (small R/a) can be corrected by 
hand for the lattice artifacts which arise since 
the lattice spatial symmetry is cubic rather than 
the continuum case which has the full rotation 
group. 

In practice, however, the b and c quarks are 
not sufficiently heavy that this static approach is 
exact. Corrections can be arranged in powers of 
1 / niQ and can be evaluated in principle using the 
heavy quark effective theory (HQET) . Examples 
of such calculations are the determination of the 
spin-orbit and spin-spin potentials between static 
quarks as well as velocity-dependent terms in the 
static potential itself. 

To explore retardation effects, one needs a 
formalism in which the heavy quarks are moving. 
One promising approach is to expand the full the- 
ory as an effective lagrangian in powers of v /c of 
the heavy quarks. This is the NRQCD scheme 
and the leading retardation effect in NRQCD 
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comes from the p. A couphng between a quark 
colour charge in motion and the gluon field. Heavy 
quark propagators are relatively easy to evaluate 
in NRQCD since the heavy quarks do not prop- 
agate backwards in time. Because there are con- 
tributions in the effective lagrangian approach of 
the form l/mqa, the continuum limit as a 
is not to be taken: instead extra terms providing 
matching with the continuum to higher powers 
in a and involving higher powers of v/c in the 
effective lagrangian are needed to increase ac- 
curacy. The coefficients of these terms should 
ideally be determined non-perturbatively but in 
practice the lowest order perturbative expressions 
(tadpole improved) are usually used. This makes 
it difficult to estimate the systematic errors in the 
NRQCD approach. Note that corrections to the 
lattice cubic symmetry to restore rotational in- 
variance will come from such higher order terms. 

Without any approximation, the lattice for- 
malism for relativistic quarks (Wilson-Dirac or 
staggered) can be applied directly to heavy quarks. 
Provided that mga << 1, this approach is quite 
straightforward. Thus only c quarks are tractable 
this way with present lattice spacings. This is a 
useful complement to the other two approaches 
which are less reliable for the case of the lighter 
c quarks. 

2. Quarkonia 

In the static approximation, a potential V{R) 
between static quarks in the fundamental rep- 
resentation of colour and at separation R can 
be extracted from the lattice. In the quenched 
approximation, this evaluation goes back to the 
early 1980's. The salient features are a behaviour 
like e/R at small R and like aR at large R. Here 
a is the string tension and e is related to the run- 
ning coupling as (indeed this is one way to deter- 
mine as from the lattice Q). The shape of the 
lattice potential (labelled S+) and the wavefunc- 
tions and energy levels are illustrated in fig. |^. 

The bb spectrum evaluated from this lattice 
potential in the Born-Oppenheimer approxima- 
tion was found not to agree precisely with ex- 
periment. One way to quantify this is that the 
energy level ratio \^Zls around 0.71 from the 
quenched lattice while it is 0.78 from exper- 



iment (for the spin averaged S and P-wave bb 
states). For a more thorough discussion of this, 
including the effect of velocity dependent terms 
in the potential, see ref 

This can be understood as a consequence of 
the quenched approximation. The value of the 
Coulomb coefficient e is expected in lowest or- 
der of perturbation theory to contain a factor 
(33 — 2Nf)~^ and so will increase as sea quarks 
(with Nf — 2, say) are included. This effect has 
been confirmed by explicit lattice calculation in- 
cluding such sea quarks 0, || . An illustration 
is shown in fig. |[ Indeed effects of including sea 
quarks seem to be comparable to those expected 
from perturbation theory although much work 
still needs to be done to include even lighter sea 
quarks in the vacuum so that the extrapolation 
to light sea quarks is under better control. The 
consequence of this increase in the depth of the 
potential at small R is that the IS level will be 
moved down in energy, resulting in an increase 
of the ratio \^Zls to bring it more into line with 
experiment. For this reason we will use differ- 
ences with the 2S energy to estimate the hybrid 
energy levels subsequently from quenched calcu- 
lations. Note that, experimentally, the IP — 15 
energy splitting is very similar for cc and bb. This 
coincidence has been used as evidence that this 
quantity is insensitive to quark masses and hence 
a good point of comparison between lattice calcu- 
lations and experiment. While this may be true 
for valence quarks, it not likely to be valid for sea 
quarks, since, as we have discussed above, the 
IS* level is especially sensitive to the sea quark 
effects. 

The other regime in which sea quarks will 
make a definite impact is in the large R region. 
It will become energetically favourable to cre- 
ate two heavy-light mesons (Qq) of energy 2mQq 
when this energy is less than V{R). This phe- 
nomenon is known as string breaking since as 
R is increased the colour flux between the static 
sources breaks with the formation of a light quark- 
antiquark pair. From the lattice mass values of 
the heavy-light mesons, this string breaking can 
be predicted to occur at around 1.2fm - an illus- 
tration is shown in fig. ^. Lattice studies of the 
potential V{R) using generahsed Wilson loops 
have not reached sufficient precision to observe 
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Figure 1: (a)Potentials and wave functions for bb hadrons in the Born-Oppenheimer approximation from 
ref 0] for quenched QCD. (b)The resulting S and P spin-averaged energy levels in the ground state potential 
(labelled E^) are compared with experiment (shown as solid lines). The energies of some of the lower lying 
hybrid excitations are also illustrated, for example in the n„ potential as will be discussed in section 3. 



this directly. It is known from studies of the ad- 
joint potential [p^ that a variational approach 
involving both the string states and the meson- 
antimeson states will be needed to obtain accu- 
rate energy estimates at these large separations 
— this is under way llJ 



14 



NRQCD calculations of quarkonia show very 
similar results to the potential approach described 
above. A comparison pO| | for spin averaged masses 
of the NRQCD result with the potential approach 
shows no significant evidence for retardation ef- 
fects. Indeed differences among NRQCD results 
[l^ arising from different lattice spacings and dif- 
ferent treatment of higher order corrections are 
of the same magnitude as their difference from 
the potential result. This is illustrated in fig. ^ 
for the quarkonium IP and 25 excitations. 

A quenched lattice study of quarkonia using 
relativistic quarks [0 also shows similar results 



to those found by the other methods described 
above. 

2.1 Quarkonia decays 

A very approximate study of some quarkonium 
decays can be made in the Born-Oppenheimer 
approximation using esentially the methods of 
atomic physics: overlaps of wave functions. In 
particular the decay to Icpton pairs will be gov- 
erned by the wave function at the origin. In 
practice the corrections to the non-relativistic ap- 
proach for decays are much larger than for the 
energy values, so this approach is rather impre- 
cise. 

Hadronic decays (such as T{AS) BB) are 
of interest because they proceed by string break- 
ing: a light quark pair is created which then re- 
sults in a pair of heavy-light mesons being pro- 
duced. This process is accessible in principle 
from lattice calculations ||l^, For example. 
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Figure 2: The potential V{R) between static quarks 
with Nf — and 2 flavours of sea quark from ref 
in units of ro ~ 0.5 fm. The lower part shows the 
differences from the common dotted curve which em- 
phasises the increase in the strength of the e/R term 
as the sea quark mass is decreased (larger k value). 



Figure 3: The potential V{R) between static quarks 
with Nf = 2 flavours of sea quark from ref Q and 
the energy corresponding to twice the ground state 
heavy-light meson mass. This shows that the string 
breaking region, where they cross, is at i? « 1.2 fm. 
No sign of string breaking is seen in this figure be- 
cause the Wilson loop operators used have a very 
small overlap with the meson-meson configuration. 



from the splitting of the energy levels caused 
by string breaking, one can estimate the decay 
rate 

3. Hybrid Mesons 

The static quark approach gives a very straight- 
forward way to explore hybrid quarkonia. These 
will be QQ states in which the gluonic contribu- 
tion is excited. The ground state of the gluonic 
degrees of freedom has been explored on the lat- 
tice, and, as expected, corresponds to a symmet- 
ric cigar-like distribution of colour flux between 
the two heavy quarks. One can then construct 
less symmetric colour distributions which would 
correspond to gluonic excitations. 

The way to organise this is to classify the 
gluonic fields according to the symmetries of the 
system. This discussion is very similar to the de- 
scription of electron wave functions in diatomic 
molecules. The symmetries are (i) rotation around 
the separation axis z with representations labelled 
by Jz (ii) CP with representations labelled by g 
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Figure 4: Excitation energies above the ground 
state for bb states from the potential approach (hor- 
izontal lines) and NRQCD (symbols) at the spatial 
lattice spacing shown in units of 2ro « 1 f m from 
for quenched lattices. The NRQCD hybrid operators 
Hi and H2 correspond to the L^*^ = 1^ and 1 
excitations respectively - these excitations are degen- 
erate in the potential approach. 

and u and (iii) GTZ. Here C interchanges Q and 
Q, P is parity and 7?. is a rotation of 180'^ about 
the mid-point around the y axis. The CTZ opera- 
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Figure 5: An illustration of the bb states in the 
static potential and in the first gluonic excitation 
from refQ using quenched lattices. Here 2ro « 1 
fm. The dotted curve illustrates the potential that 
would be needed to reproduce more closely the ob- 
served spectrum of 15", IP and 2S states. 

tion is only relevant to classify states with = 0. 
The convention is to label states of Jz — 0, 1, 2 
by E, n, A respectively. 

In lattice studies the rotation around the sep- 
aration axis is replaced by a four-fold discrete 
symmetry and states are labelled by represen- 
tations of the discrete group D^h- The ground 
state configuration of the colour flux is then S+ 
{Aig on the lattice). The exploration of the en- 
ergy levels of other representations has a long 
history in lattice studies Q|. The first excited 
state is found to be the !!„ on a lattice) - see 
fig. H for an illustration. This can be visualised 
as the symmetry of a string bowed out in the 
X direction minus the same deflection in the —x 
direction (plus another component of the two- 
dimensional representation with the transverse 
direction x replaced by y), corresponding to flux 
states from a lattice operator which is the differ- 
ence of U-shaped paths from quark to antiquark 
of the form □ — U. 

Recent lattice studies have used an asym- 
metric space/time spacing which enables excited 
states to be determined in a well controlled way. 
Results are shown in fig. |6| for a large variety 
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Figure 6: The potential energies, extrapolated to 
the continuum limit, for different gluonic excitations 
from refjl^] for quenched lattices. Here 2ro ~ 1 fm. 

of gluonic excitations. These results conflrm the 
finding that the IIu excitation is the lowest lying 
and hence of most relevance to spectroscopy. 

From the potential corresponding to these 
excited gluonic states, one can determine the spec- 
trum of hybrid quarkonia using the Schrodinger 
equation in the Born-Oppenheimer approxima- 
tion. This approximation will be good if the 
heavy quarks move very little in the time it takes 
for the potential between them to become es- 
tablished. More quantitatively, we require that 
the potential energy of gluonic excitation is much 
larger than the typical energy of orbital or radial 
excitation. This is indeed the case especially 
for b quarks. Another nice feature of this ap- 
proach is that the self energy of the static sources 
cancels in the energy difference between this hy- 
brid state and the QQ states. Thus the lattice 
approach gives directly the excitation energy of 
each gluonic excitation. 

The n„ symmetry state corresponds to exci- 
tations of the gluonic field in quarkonium called 
magnetic (with L^^ — 1^ ) and pseudo-electric 
(with I ^) in contrast to the usual P-wave or- 
bital excitation which has L^'-' = 1 . Thus 
we expect different quantum number assignments 
from those of the gluonic ground state. Indeed 
combining with the heavy quark spins, we get a 
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degenerate set of 8 states with J^^ = 1 , +, 
2"+ and 1++, 0+", 2+" respectively. 

Note that of these, J^^ = 1"+, and 2+" 

are spin-exotic and hence will not mix with QQ 
states. They thus form a very attractive goal 
for experimental searches for hybrid mesons. Il- 
lustrations of the spectrum of such spin-exotic 
hybrid mesons are given in fig. and |. 

The eightfold degeneracy of the static ap- 
proach will be broken by various corrections. As 
an example, one of the eight degenerate hybrid 
states is a pseudoscalar with the heavy quarks in 
a spin triplet. This has the same overall quantum 
numbers as the S-wave QQ state {rjb) which, how- 
ever, has the heavy quarks in a spin singlet. So 
any mixing between these states must be medi- 
ated by spin dependent interactions. These spin 
dependent interactions will be smaller for heavier 
quarks. It is of interest to establish the strength 
of these effects for b and c quarks. Another topic 
of interest is the splitting between the spin exotic 
hybrids which will come from the different ener- 
gies of the magnetic and pseudo-electric giuonic 
excitations. 

One way to study this is using the NRQCD 
approach which enables the L^'~'' = 1^ and 
1 excitations to be separated in a spin aver- 
aged approach. Lattice results indicate no 
statistically significant splitting (see fig. ||) al- 
though the 1^ excitation does lie a little lighter. 
This would imply, after adding in heavy quark 
spin, that the J^*-^ = 1 ^ hybrid was the light- 
est spin exotic. In principle the NRQCD ap- 
proach, by adding spin-dependent terms in the 
Lagrangian, can address the full sphtting of the 
8 levels. However, as we shall also discuss in con- 
nection with propagating quark approaches, the 
mixing of non spin exotic states with QQ may 
confuse this situation. Including spin-dependent 
terms in a NRQCD study of hybrids does give |^ 
a relatively large spin splitting among the triplet 
states. Unfortunately this study has only consid- 
ered magnetic giuonic excitations so cannot ad- 
dress the splitting between spin exotic hybrids. 

Confirmation of the ordering of the spin ex- 
otic states also comes from lattice studies with 
propagating quarks jl^, |l^ which are able to 
measure masses for all 8 states. We discuss this 
evidence in more detail below - see also fig. 0. 



Within the quenched approximation, the lat- 
tice evidence for bb quarks points to a lightest 
hybrid spin exotic with J^*^ = 1 ^ at an en- 
ergy given by (mjy — TO2s)''o =1.8 (static poten- 
tial [|); 1.9 (static potential @, NRQCD [|); 

2.0 (NRQCD These results can be sum- 
marised as 

(mn - m2s)ro = 1.9 ± 0.1 

Here ro is defined implicitly by the static force as 
r'^F{r) — 1.65 at r = tq and is a well measured 
quantity on a lattice derived from the static po- 
tential V{R) at r « 0.5 fm. Within the quenched 
approximation, where different experimental ob- 
servables differ by of order 10%, the overall scale 
is uncertain but we choose r^^ = 390 MeV ± 
10%. Using the experimental mass of the T(2S'), 
this implies that the lightest spin exotic hybrid is 
at iriH — 10.76(7) GeV. Above this energy there 
will be many more hybrid states, many of which 
will be spin exotic. 

Some preliminary lattice studies have been 
made including sea quarks. As yet only sea quark 
masses down to the strange quark mass have 
been explored and hence the extrapolation to re- 
alistic sea quark masses is not yet well estab- 
lished. Indeed, pushing to lower sea quark masses 
is the main remaining challenge in lattice gauge 
theory. The light propagating quark case has 
been explored ||l^ but no significant differences 
are found from using quenched vacua. The ex- 
cited giuonic static potential has also been de- 
termined including sea quarks {Nf = 2 flavours) 
and no significant difference is seen iQ. Thus 
the quenched estimates given above are not su- 
perseded. Note, however, that hybrid states can 
mix with QQqq states - for instance with their 
decay products as we shall discuss in the next 
section. This mixing is, in principle, enabled in 
a lattice study with sea quarks. 

3.1 Light quark hybrid mesons 

Here we focus on lattice results for hybrid mesons 
made from light quarks using fully relativistic 
propagating quarks. There will be no mixing 
with qq mesons for spin-exotic hybrid mesons and 
these are of special interest. The first study of 
this area was by the UKQCD Collaboration |l^ 
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Figure 7: The masses in GeV of states of J^'^ built 
from hybrid operators with strange quarks, spin- 
exotic (*) and non-exotic (squares). The dot-dashed 
hnes are the mass values found for ss operators. 
Quenched lattice results from ref rl6l . 



who used operators motivated by the heavy quark 
studies referred to above. Using non-local op- 
erators, they studied all 8 J^^ values coming 
from L^'^ = 1"' and 1 ^ excitations. The re- 
sulting mass spectrum is shown in fig. ^ where 
the J^'~'' = 1 ^ state is seen to be the lightest 
spin-exotic state with a statistical significance of 
1 standard deviation. The statistical error on 
the mass of this lightest spin-exotic meson is 7% 
but, to take account of systematic errors from 
the lattice determination, a mass of 2000(200) 
MeV is quoted for this hybrid meson with ss light 
quarks. Although not directly measured, the cor- 
responding light quark hybrid meson would be 
expected to be around 120 MeV lighter. 

One feature clearly seen in fig. ^ is that non 
spin-exotic mesons created by hybrid meson op- 
erators have masses which are very similar to 
those found when the states are created by qq op- 
erators. This suggests that there is quite strong 
coupling between hybrid and qq mesons even in 
the quenched approximation. This would imply 
that the 7r(1800) is unlikely to be a pure hybrid, 
for example. 

A second lattice group has also evaluated hy- 
brid meson spectra with propagating quarks from 
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Figure 8: The extrapolation in sea quark mass (in 
lattice units) from ref[^^ for the 1 hybrid me- 
son. The dotted vertical line corresponds to light sea 
quarks. 

quenched lattices. They obtain Q masses of the 
1 state with statistical and various systematic 
errors of 1970(90)(300) MeV, 2170(80)(100)(100) 
MeV and 4390(80) (200) MeV for nn, ss and cc 
quarks respectively. For the 0'' spin-exotic state 
they have a noisier signal but evidence that it is 
heavier. They also explore mixing matrix ele- 
ments between spin-exotic hybrid states and 4 
quark operators. 

Recently a first attempt has been made 
to determine the hybrid meson spectrum using 
full QCD. The sea quarks used have several dif- 
ferent masses and an extrapolation is made to 
the limit of physical sea quark masses, yielding 
a mass of 1-9(2) GeV for the lightest spin-exotic 
hybrid meson, which they again find to be the 
1 '". In principle this calculation should take ac- 
count of sea quark effects such as the mixing be- 
tween such a hybrid meson and qqqq states such 
as rjTT. As illustrated in fig. ||, the calculations are 
performed for quite heavy sea quarks (the light- 
est being approximately the strange quark mass) 
and then a linear extrapolation is made. It is 
quite possible, however, that such mixing effects 
turn on non-linearly as the sea quark masses are 
reduced. The systematic error from this possi- 
bility is difficult to quantify. 

The three independent lattice calculations of 
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the hght hybrid spectrum are in good agreement 
with each other. They imply that the natural en- 
ergy range for spin-exotic hybrid mesons is around 
1.9 GeV. The J^'^ = 1"+ state is found to be 
lightest. It is not easy to reconcile these lattice 
results with experimental indications for res- 
onances at 1.4 GeV and 1.6 GeV, especially the 
lower mass value. Mixing with qqqq states such 
as TyTT is not included for realistic quark masses in 
the lattice calculations. This can be interpreted, 
dependent on one's viewpoint, as either that the 
lattice calculations are incomplete or as an indi- 
cation that the experimental states may have an 
important meson-meson component in them. 

3.2 Hybrid meson decays 

One clear feature of heavy quark hybrid mesons 
is that they have very extended wavefunctions 
since the potential that binds them is relatively 
flat. This has implications for their production 
and decay. For instance, any vector state will 
only be weakly produced in e^e^ collisions be- 
cause the wave function at the origin will be 
small. 

Given our mass estimates above, the open 
channels for decay of a J^'~"' = 1 ^ hybrid in- 
clude BB,BB*,r]br],r]br]',T{lS)uj and T{lS)<p. 
Selection rules have been proposed for hybrid de- 
cays, for example |2l[ that H X + Y if X and 
Y have the same non-relativistic structure and 
each has L = 0. This would rule out BB and 
BB* and the analogous cases for charm quarks. 

This selection rule can be addressed directly 
from the static quark approach. The symme- 
tries in this case of rotations about the separation 
axis, etc have to be preserved in the strong decay. 
From the initial state with the giuonic field in a 
given symmetry representation, the qq pair must 
be produced in the decay in such a way that the 
combined symmetry of the quark pair and the 
final giuonic distribution matches the initial rep- 
resentation. 

For the ground state of the giuonic excita- 
tion (non- hybrid) we have Jj, = and even CP. 
Thus, for this state to decay to {Qq){Qq) with 
each heavy-light meson having L = 0, the final 
giuonic distribution is also symmetric (actually 
it is essentially two spherical blobs around each 
static source binding the heavy light mesons). 



Then any qq pair production has to respect this 
symmetry and have = and even CP. Since 
there is no orbital angular momentum, the CP 
condition then requires Sqq = 1, a triplet state. 
This is just a derivation of what is called the 
^Po model of decays: the light quark-antiquark 
is produced in a triplet state. This spin assign- 
ment can be tested by the ratio of BB, BB* and 
B*B* decays. 

For the J^^ = 1 ^ hybrid we have a giuonic 
field with = 1 and odd CP. For the case of 
decay to a {Qq){Qq) with each heavy- light me- 
son having L = 0, this would imply that the 
qq would have to be produced with = 1 and 
odd CP. This is not possible since the triplet 
state would have even CP while the singlet state 
cannot have = 1. This is then equivalent to 
the selection rule described above. There will 
presumably be small corrections to this selection 
rule coming from retardation effects. Decay to 
{Qq){Qq) with one heavy- light meson having a 
non-zero orbital excitation is allowed from sym- 
metry but is not allowed energetically with con- 
ventional mass assignments for the P-wave ex- 
cited B meson multiplet. 

Decays to {QQ){qq) are also possible since 
there is enough excitation energy to create a light 
quark meson. This meson must be created in a 
flavour singlet state and the lightest candidates 
are 77 and uj. In a lattice context, this produc- 
tion is via a disconnected quark loop with s, u 
and d quark contributions of similar strength. 
The flavour singlet mixture of rj and 77' (mainly 
r/) and the singlet mixture of the vector mesons 
(which includes a substantial w component) are 
expected to be coupled most strongly. 

So allowed decays are rjbr], ij^r]' , T{1S)uj and 
T{lS)(j). Here the light meson must have Jz = 1 
and together with the CP constraint, this implies 
that the light meson must be in a P-wave with 
respect to the heavy quark meson. 

S'-wave decays to rjbfi and T{lS)h are also 
allowed although there may be insufficient phase 
space. (Here the /i and h are J^'~'' = 1++ and 
1 ^ flavour singlet mesons). 

As for the case of quarkonium decays, it is 
possible in principle to explore on the lattice some 
aspects of these decays. One can study matrix 
elements between ground states which are degen- 
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erate in energy such as the 1 ^ hybrid and the 
ribT] final state where the light quark mass is ad- 
justed so that there is equal energy in both sys- 
tems. This and similar lattice studies will enable 
some further guidance to be given for experimen- 
tal searches for hybrid mesons. 

4. Summary and Outlook 

One of the advantages of lattice studies is that, 
by varying the parameters such as quark masses, 
they can serve as very useful data to develop 
phenomenological models. One example is that 
the excitation spectrum of the potential between 
static quarks can be used to test QCD string ex- 
citation models. This has been much discussed - 
for a review see ref . 

At present, lattice studies are restricted to 
sea quark masses no lighter than strange quarks. 



ods and lie around ^pj^fg w 3.0 but the system- 
atic errors are quite large. 

The situation is much better controlled for bb 
hybrids. We expect the lightest bb hybrid to have 
J^^ = 1-+ and mass 10.75 ± 0.10 GeV. It will 
be difficult to isolate such states experimentally - 
but well worth the effort. The most likely decay 
modes of this hybrid meson are to a bb ground 
state meson {rjb or T(IS')) with the emission of 
a flavour singlet light quark meson (77, 77', uj or 
(/)). Future lattice calculations should be able to 
study these and other decays. 
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